function plot_running_modulation(proj_meta,cur_site,cell)
runM=proj_meta(cur_site).rd(1,1).velM_smoothed>0.01;
if cell>0
    act=proj_meta(cur_site).rd(1,1).act(cell,:);
    figure('color','k');
    set(gcf,'pos',[17         710        3804         273])
    area([1:length(runM)],5*[runM],'facecolor',[0 0.8 1],'linestyle','none')
    hold on
    plot(act,'k')
    box off
    axis tight
    ylim([0.5 4])
else
    act=proj_meta(cur_site).rd(1,1).act;
    figure('color','k');
    set(gcf,'pos',[17         710        3804         273])
    area([1:length(runM)],5*[runM],'facecolor',[0 0.8 1],'linestyle','none')
    hold on
    plot(act','k')
    box off
    axis tight
    ylim([0.5 4])
end